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^ ■ Abstract 



We obtain analytic, nonperturbative, approximate solutions of Yukawa theory in 
the one-fermion sector using light-front quantization. The theory is regulated in 
the ultraviolet by the introduction of heavy Pauli- Villars scalar and fermion fields, 
each with negative norm. In order to obtain a directly soluble problem, fermion-pair 
creation and annihilation are neglected, and the number of bosonic constituents is 
limited to one of either type. We discuss some of the features of the wave function 
of the eigensolution, including its endpoint behavior and spin and orbital angular 
momentum content. The limit of infinite Pauli- Villars mass receives special scrutiny. 
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1 Introduction 



It is remarkable that the first analysis of renormalization in quantum field 
theory was performed in the context of a nonperturbative problem - the QED 
calculation of the Lamb Shift in hydrogenic atoms by Bethe in 1948 [1]. In 
effect, one calculates the energy shift to order a and all orders in the atomic 
binding from the difference [2]: 

AE n =< ^\T,(Za) - bm\i\) n > , (1) 

where E(Za) is the order-a self energy of the electron evaluated in the back- 
ground of the nuclear Coulomb potential and 8m = S(Z = 0) is the free- 
electron mass shift. The —5m = m — m counterterm arises from the order-a 
contribution to the physical mass. As shown by Yennie and Erickson [2], the 
atomic energy shift AE n can be evaluated in terms of the expectation values 
of gauge invariant operators which are functionals of the field strength F^ v of 
the background field. The dominant term in the shift of the nS levels is of 
order 

AE n ~ —{Zaf\Yv(Zafm. (2) 

7T 

The logarithmic dependence indicates the inherently nonperturbative nature 
of the renormalization problem. The atomic-physics analysis of renormaliza- 
tion is appropriate to cases such as an infinitely heavy nucleus where the 
binding interaction is effectively static. In principle, one can take into account 
dynamical effects from retarded interactions and finite mass sources by using 
effective interaction methods [3]. 

Relativistic systems in quantum field theory can be analyzed by quantization 
on the light front [4] in a manner resembling nonrelativistic theory. The eigen- 
states of the invariant light-front Hamiltonian H LF = P~P + — Pi satisfy the 
LF Heisenberg equation 

H hF \^j n >= M*\ij n > . (3) 

The eigenvalue problem can be rewritten in matrix form by introducing a 
free Fock basis. The matrix equation becomes discrete in momentum space by 
introducing periodic boundary conditions, as in the discrete light-cone quan- 
tization (DLCQ) method [5,6]. The theory can be rendered ultraviolet finite 
by introducing Pauli-Villars (PV) ghost fields [7,8,9,10,11,12,13]. We have 
recently demonstrated the viability of this type of PV regularization by ap- 
plying it to (3+l)-dimensional Yukawa theory of spin-half fermions and scalar 
bosons and have obtained nonperturbative DLCQ solutions of this theory in 
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low particle number sectors [10]. We have also been able to solve such theories 
analytically in the limit of exact degeneracy of the negative and positive norm 
states [11]. Infrared divergences do not appear in neutral bound states such 
as color-singlet hadrons in QCD, since the gauge-particle interactions cancel 
at long wavelength. 

The interaction terms in H^p change particle number and contain the quantum 
effects associated with dressing the constituents analogous to T,{Za). The 
counterterms arising from the difference of the physical and bare mass of the 
massive constituents give a subtraction term [14] analogous to the subtraction 
term in Eq. (1), which is 

_ §Tn? 

SM 2 n = -J2<^n\^ L \^n>, (4) 
q x q 



where — = ^ is evaluated within the sum over Fock states. 

In practice, carrying out the above analysis is complicated in nonperturbative 
relativistic problems by the necessity that the Fock state has to be truncated. 
However, we do not perform sector-dependent renormalization [15], and in- 
stead retain the bare parameters of the original Lagrangian, the fermion mass 
and the coupling, as the only parameters to be adjusted. The shift in the 
fermion mass is the only counterterm. Contributions to the mass shift from 
instantaneous fermion interactions do not appear if one includes a negative- 
metric fermion as part of the Pauli-Villars regularization [11], since these 
instantaneous terms are mass independent. 

In this paper we shall carry out the above nonperturbative mass renormal- 
ization program in the context of the LF Hamiltonian formulation of Yukawa 
theory. The interactions of a massive fermion and light scalar boson generate 
an effective fermion bound state, for which we are able to obtain approximate, 
but analytic, nonperturbative solutions. A set of negative-norm heavy scalar 
and fermion Pauli-Villars fields is used to regulate the ultraviolet divergences 
while preserving the chiral symmetry of the perturbative mass shift. In order 
to obtain a directly soluble problem, fermion-pair creation and annihilation 
are neglected, and the number of bosonic constituents is limited to one of 
either type. We shall discuss some of the features of the wave function of the 
eigensolution, including its endpoint behavior and spin and orbital angular 
momentum content. The dependence of the renormalized theory on the mass 
of the Pauli-Villars fields has some unexpected features which we discuss in 
detail. 

In principle, the light-front wave functions of QCD can be computed directly 
by the diagonalization of the light-front Hamiltonian. The DLCQ method [5,6] 
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provides a discretization scheme which transforms the eigenvalue problem 
of QCD into the problem of diagonalizing very large sparse matrices. Al- 
though the Fock space is truncated, the DLCQ method retains the essential 
Lorentz symmetries of the theory including boost independence. DLCQ has 
also provided an important tool for analyzing string and higher dimension 
theories [16]. The DLCQ method has been successfully applied to QCD and 
other gauge theories in one space and one time dimensions [6]. There have 
also been applications of DLCQ to (3+l)-dimensional non-gauge theories [6]. 

Application to QCD 3+1 is computationally intensive [17], because of the large 
numbers of degrees of freedom; however, the solution of the bound-state eigen- 
value problem corresponding to the hadronic spectrum would be a very im- 
portant step. Given the projection of the eigensolutions on the light-front Fock 
basis, one can compute observables such as form factors and transition am- 
plitudes [18,19], the underlying features of deep inelastic scattering structure 
functions, the distribution amplitudes which control leading twist contribu- 
tions to hard exclusive processes [20], and the skewed parton distributions 
which can be measured in deeply virtual Compton scattering [21,22]. First- 
principle computations of exclusive decay amplitudes of heavy hadrons, such 
as the D and B mesons [23,24], require knowledge of heavy and light hadron 
wave functions in order to extract the phases and other parameters of the 
electroweak theory. Light-front techniques can also be applied to traditional 
nuclear physics [25]. The light-front representation is boost-independent and 
provides the nonperturbative input and matrix elements required for such 
analyses. 

An important feature of light-front Hamiltonian is the simplicity of spin and 
angular momentum projections: the sum rule for the angular momentum of 
the eigensolution J z = Yn=i Sf + J2?=i M holds Fock state by Fock state. 
Here the sum is over the spin projections S z of the constituents in the n- 
particle Fock state. There are only n — 1 contributions to the internal orbital 
angular momentum. The spin projections also provide a convenient way to 
classify independent contributions to the wave functions. Our truncation to 
two partons limits L\ to the values and ±1. 

In Sec. 2 we discuss the Yukawa Hamiltonian and its regular izat ion and renor- 
malization, as well as the Hamiltonian eigenvalue problem, which we solve, 
given the truncation to two particles. Since the methods we use are analytic, 
we are able to find two-parton solutions in the continuum theory without 
DLCQ or other discretization. We discuss the nature of the solutions in two 
limits, one in Sec. 3.1 where the PV masses are equal and another in Sec. 3.2 
where the PV boson mass approaches infinity more slowly than the PV fermion 
mass. We are particularly interested in the chiral properties, the large trans- 
verse momentum fall-off, and the end-point behavior of the eigensolutions in 
the light-cone variables Xi = k + / P + of the constituents. In Sec. 4 we argue 
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that in calculations where the representation space is truncated, it is nec- 
essary to keep the values of the PV masses finite even in cases where it is 
computationally possible to take the limit of infinite PV masses. It is possible 
that some of the effects we see may be related to the triviality of Yukawa 
theory [26], but similar considerations probably apply to asymptotically free 
theories. Section 5 contains our conclusions. 

Our calculations are somewhat similar to those of Bylev, Glazek, and Przes- 
zowski [27], except that they did not use a covariant regulation procedure; it 
is the effect of the covariant regulator that will form the focus of our discus- 
sion here. Similar work in a purely scalar theory has been done by Bernard et 
al. [28]. For a more formal treatment of dressed constituents, see [29]. 

The notation that we use for light-cone coordinates is 

x ± = x° + x 3 , x± = (x 1 ,x 2 ). (5) 

The time coordinate is x + , and the dot product of two four- vectors is 

p ■ x = -(p + x~ + p~x + ) — p± ■ x± . (6) 



The momentum component conjugate to x is p + , and the light-cone energy 
is p~ . Light-cone three- vectors are identified by underscores, such as 

p = (p + ,p ± ). (7) 



For additional details, see Appendix A of Ref . [8] or the review [6] . 



2 Computational Framework 



Taking the physical fermionic and bosonic fields to be ipi and 0i, respectively, 
and the PV (negative-metric) fields to be ip2 and 02, the Yukawa action be- 
comes 



d A x 



(8) 
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where the scalar three-point interaction is expressed in terms of zero-norm 
fields 



i> = (V>i + ih) , (j>= (0i + <h) • (9) 

For simplicity, we will not consider a 4 term; with pair creation removed, it 
will not be required. 

If the bare mass of the fermion is zero, the bare action (the same as above but 
with the PV fields set equal to zero) possesses a discrete chiral symmetry - 
invariance under the transformations ip — > 75^, — > —0. One consequence of 
this symmetry is that the physical mass is also zero. The full action, including 
the PV fields breaks this symmetry explicitly. There are two possible versions 
of this breaking: if we require that the PV Fermi field is not transformed, 
the cross term in the interaction breaks the symmetry; if we require the PV 
Fermi field to transform in the same way as the physical field, the PV mass 
term breaks the chiral symmetry. Thus an interesting question, which we will 
examine below, is whether chiral symmetry will be restored in the limit of 
large PV masses where the unphysical states decouple. It will be a point of 
interest to see if the symmetry is restored, at least in the sense that the bare 
mass and the physical mass are proportional to each other. 

The corresponding light-cone Hamiltonian, except for the addition of the PV 
fields, has been given by McCartor and Robertson [30]. Here we include the PV 
fields but neglect pair terms and any other terms which involve anti-fermions. 
The resulting Hamiltonian is 



l.S r 



(10) 



+ E/d9^^(- 1 ) i+1 «5(9K(9) 
3 q 

+ J2 J d P d l { [ V -2 S (P, 9) + v 2s(p + g, qj\ bl s (p)al(q)bi_ s (p + q) 



i,j,k,s 



+ [Uj(p,q) + Ui(p + q,qj\ b] :S (p)al(q)b iyS (p + q) + h.c.} , 



where 

Uj(p,q) = 



V167T 3 P+Vq 1 



V^p + V¥' 



(ii) 



mi is the mass of the bare fermion, jj, = [i\ is the physical boson mass, [12 
and m 2 are the masses of the PV boson and fermion, respectively, and e*2s = 
— 4=(2s, «). The V interaction introduces one unit of relative orbital angular 
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momentum projection L z which is compensated by the change in fermion spin 
projection S z to conserve J z . [31]. The nonzero commutators are 



[ai(9),o}(g')]=(-l) i M(9-2'), 

{bi,s(p), b{ s ,{p)} = (-mAs'Sip ~ g) ■ (12) 

The Fock-state expansion for a spin-1/2 fermion eigenstate of the Hamiltonian 
is 

OO ntot n _ fctot „ "tot fctot 

**= E Il/^EIl/^^-E^-Ea*) (is) 

Tii,n2,fci,fe2=0n=l fc=l n fc 

1 "tot fctot 

*fc fcj) (p n ;g fc ) 7= — 1 n *U(ej n 4(2 fc )io) , 

yIl»Wi!lljKj! n fc 

where ni is the number of bare fermions, n 2 the number of PV fermions, k\ 
the number of physical bosons, and h 2 the number of PV bosons. The total 
number of fermions is n tot = rt\ + n 2 , and the total number of bosons is 
fctot — fci + fc2- The n-th constituent fermion is of type i n , and the fc-th boson 
is of type jk- This Fock state expansion will be used to solve the eigenvalue 
problem (P + P~ -P^)$ CT = M 2 $ CT . The normalization of the eigenstate is 

&*-&* = 5(P!-P). (14) 
Our first approximation will be to truncate the expansion to two particles: 

= E *4f(£)|o> + E / dij ijs (i_)bl(p - iWmo) (is) 



and reduce the eigenvalue problem by projecting onto Fock sectors. Without 
loss of generality, we consider only the J z = +1/2 case. The resulting coupled 
equations determine the wave functions fcj± to be 



fi j+ (L) 



P+ 

1V1 1-1+ /P+ 1+/P+ 



X 



Q2(-i) k+1 z k )Ui(p + Y,(-i) k+1 z k u k (p,i) 
ik k 

J^2 _ m i +i ± M 3 +' ± 



(16) 



i-i+/pi 
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corresponding respectively to Fock components where the fermion constituent 
is aligned or anti- aligned with the total spin J z . The nonperturbative physics 
is thus contained in the determination of Zi and M 2 . When the two-body wave 
functions / are eliminated from the one-fermion projections, we obtain 1 

(M 2 - m 2 )zi = g 2 fij(zi - z 2 )J + g 2 m i (zim 1 - z 2 m 2 )I 

+g 2 /i 1 [(z 1 - z 2 )mi + z 1 m 1 - z 2 m 2 \I 1 , (17) 

with 



J 



dydl\ 


E 

jk 




(-iy +k 


(m j /m) n 


16tt 2 


M 2 - 


m )+i\ fat* 


y(i-y)"' 




i-y y 


dydl\ 


E 

jk 




(_!)*+* 


(m 2 + Z 2 )//i 2 


167T 2 


M 2 - 




y(i -y) 2 




1-2/ y 



i 16tt 2 4^ ^2 y(l -y) 2 ' { ' 



m 2 +l 2 a 2 +l 2 

These integrals are not independent; a change of variable tow — j' x + y ± 
and an interchange of integration order can be used to show that 

M 2 

J = -5-Jo • (20) 
pi 



We solve the i = 2 case of Eq. (17) for ( = Z2/21, to obtain 



g 2 /x 2 J + g 2 fii(m 1 + m 2 )/i + g 2 m 1 m 2 I { 



~~ M 2 — m 2 + 5f 2 /i 2 J + 2g 2 fi 1 m 2 I 1 + 2 m| J ' 1 j 



From the remaining i = 1 case we solve for g . This yields 



., (MT ffll )(MTm 2 ; 



, = v : : ti_ / 22 n 

y (m 2 -m 1 )( f i 1 I 1 ±MI ) 1 ; 



There are two possible solutions for g 2 since the remaining equation is quadratic 
in g 2 . Substitution into (21) and use of (20) reduces ( to the remarkably simple 
form 

C = MTmi , (23) 



independent of I and I\. 

1 Note that the projection onto the opposite spin is automatically zero because the 
integrand is linear in l±. 
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We pause here to remark that the wave function we have obtained in the 
nonperturbative calculation is very similar to the one which we would obtain 
using first-order perturbation theory to perturb about the state of one bare, 
physical fermion. The only differences are that in perturbation theory z 2 = 
and M = mi. From (23) we see that as m 2 — > oo, ( = z 2 /zi will be small as 
long as mi << m 2 . Since this last requirement is necessary if we are to expect 
to restore at least approximate unitarity in the limit of large PV masses, we 
will insist on it. The only significant difference between our nonperturbative 
calculation and first-order perturbation theory is that in perturbation theory 
M = mi while in the nonperturbative calculation M is determined by (22). 2 
The form of the wave function in terms of the parameters is exactly the same 
in the perturbative and nonperturbative calculations; only the parameters are 
different. 

In the presence of the negatively normed constituents, we define the "physical 
wave functions" as the coefficients of Fock states containing only positive-norm 
particles. This can be done without ambiguity by requiring that all Fock states 
be expressed in terms of the positive-norm creation operators b\ s and a\ and 
the zero-norm combinations b\ = b\ s + b 2s and eft = a[ + a\. Because b\ is 
null, a fermion created by b\ is annihilated by the generalized electromagnetic 
current + ip 2 )^{jpi + V^) appropriate to this PV-regulated theory; thus 
the null fermions do not contribute to current matrix elements and should 
not make a physical contribution to a state. By analogy, eft is also deemed 
to create unphysical contributions. The procedure, then, is to express the 
wave function in terms the operators b\ s , b\, a\ and aft acting on the vacuum. 
Any term containing a b\ or an oft is then discarded when constructing the 
physical state. This procedure is a non-orthogonal projection onto the physical 
subspace. 

After application of this procedure to our case, the physical state with spin 
J z = +1/2 in the two-particle truncation is 



<S> { +ly S = (Zl-Z 2 )b[ + (P)\0) 

+ £ I dq(-l) i+ i f ijs (q)b\ s (P - q)a\(q)\0) . (24) 

i,j,s 

The normalization condition (14) fixes Z\. In addition to fixing the physi- 
cal mass, one additional renormalization condition is needed. In previous pa- 
pers [8,9,10] we have specified a value for the expectation value in the state 

(:0 2 (O):> = $t;0 2 (O):$ CT . (25) 

2 In practice we will fix M as a renormalization condition and use (22) to restrict 
the behavior of g and mi as functions of m 2 . 
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For some of the solutions given below this quantity diverges even after renor- 
malization, so this is not a suitable condition. In the rest of the paper it is not 
necessary to specify the final renormalization condition. In place of specifying 
the final normalization condition, we will examine features of the solution. We 
will look for cases where the structure functions are finite and nonzero. 



The normalization of $+phy S an d the definition of (:0 2 (O):) reduce to 



i = (^-^) 2 + E / dl - 



E(-!) l+ v^a) 



■ij 



! (0):)=E/^ 



1+/P+ 



E(-!) l+ %*a) 



(26) 
(27) 



These relations can be written more explicitly in terms of the following inte- 
grals: 



rlti ~ f rln 9 



where 



(28) 



it) 



(y) = J dil E 



h(y) 



k M 2 — 



y 



lA 



dl 2 ± 
( 



Uli 



X 



E 



(-iy +k 

M 2_^l+i_ I y(l - y) 



1-2/ 



j (y) = f dl 2 ± IE ( ^IT^ 2 2 

V i-?/ y 



\ 



(29) 
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3l {y) = /< E 



(-iy +k 



\ 



i-?/ y 



For the normalization and for (:0 2 (O):), we then obtain 



-, = (i-(Y[i + gVo + Ji)] 



(:0 2 (O):) = ^ 2 {(1-C) 2 [4 + J(] 



(mi - Cm 2 ) 2 



^0 + 2(1 - C) 



4 ■ 



(30) 



(31) 



The boson structure functions are given by 



d 2 l 



(32) 



In terms of integrals already defined in (29) we obtain 



fB+(y) = 



16tt 2 



+ 2(i- C ) ( ""- Cm2 W 
/B-(y) = f^(i -C) 2 ji(y)- 



(33) 
(34) 



As an alternative renormalization condition one could use the radius R of the 
dressed- fermion state, as defined by the slope of the Dirac form factor F ± . 
These quantities are related by the standard expression R = ^—6F{(0). The 

slope can be computed from the eigenfunction $+p hys as 



I? 
6 



F{(0) = zl U 



16tt 2 



i 

£(-iy' +/ E(-iy + V a(a-l)day 3 dy 

i'j' ij o 



(35) 



x 



{ (1 - Q 2 mim v + (mi + m;/)(l - C)("ii - C^)(l - y) 



+ (m 1 -Cm 2 ) 2 (l-y) 2 l /{Dtff 
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+ 2(1-07^!?'}, 



with Dff = a(ym 2 + (1 - y)/i 2 ,) + (1 - a)(ym 2 + (1 - y)/i 2 ) - y(l - y)M 2 . 
Similarly we can extract the axial coupling [32] 



9A = zl(l- C) 2 + zf^- 2 E(-!) l+J / dctydy (36) 

i'j' ij o 

x | (1 - () 2 mim v + (m, + m v )(\ - Q{ m i ~ C m 2)(l - y) 

+ (m 1 -(m 2 ) 2 (l-y) 2 } 

+ (l-C) 2 log[2Di/]} 
and the anomalous magnetic moment k = F 2 (0) of the dressed fermion 



K = 2M^ 2 (l-C) T f^E(-ir' +/ E(-l) l+ V day 2 dy (37) 



1 j u 



x [(1 - C) (am f + (1 - a)mi) + (mi - C^)(l - y)\ /D. 



The result for the anomalous moment is confirmed by comparison with Eq. (51) 
of Ref. [31]. If the fermion x in Ref. [31] is written as 1 — y and the M in the 
numerator is replaced by mi, as per the discussion after Eq. (46), the two 
results agree, once we drop the sum over PV particles. Note that only the 
two-particle Fock state contributes since the anomalous moment requires a 
change in L z without a change in particle number. 

As an example of how R might be used as a renormalization condition, we 
compute R and g 2 for a series of mi values, with m 2 — fJ>2 — 10/ii and M 
fixed at /ii. 3 We use the lower signs in Eqs. (22) and (23). The results are 
plotted in Fig. 1. The figures show that for a chosen value of R and M one can 
obtain values for the bare parameters g and mi, with the only ambiguity being 
between weak and strong coupling. The axial coupling is essentially constant 
over the given range, at a value very near unity. The anomalous moment is 
plotted in Fig. 2; here there is again the double- valued structure of the radius 
R. That the anomalous moment vanishes as MR — > is consistent with the 
Drell-Hearn-Gerasimov sum rule, as discussed in [32]. 



The boson is not dressed due to the fact that we have eliminated pair production. 
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Fig. 1. Plots of (a) the dressed-fermion radius R and (b) the bare coupling squared 
g 2 as functions of the ratio of the bare fermion mass mi to the physical boson mass 
fii. The PV masses are fixed at mi = fi2 = 10/xi, and the dressed-fermion mass at 
M = m. 




Fig. 2. The anomalous moment k of the dressed fermion, multiplied by 10 2 , as a func- 
tion of its radius R, scaled by its mass M. For this particular plot the dressed-fermion 
mass is set equal to the physical boson mass fii, and the PV masses are fixed at 
1^2 = A*2 = 10/ii. The solid lines are drawn to connect points on the same branch. 
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3 Limits for large Pauli-Villars Masses 



3.1 Equal Pauli-Villars Masses 

Although the nonperturbative problem has been reduced to a single nonlinear 
equation, and although all the integrals involved in that equation can be done 
in closed form, the resulting expression is very long and complex and, worse 
yet, is a function of many variables. To gain some control over the total space 
in which we will look for solutions, we will fix the ratio of the two PV masses. 
A natural choice seems to be 7712 = fJ>2, especially if we choose M, the physical 
fermion mass, to be equal to //1, the physical boson mass. 

When the PV masses m 2 and /x 2 are equal, and we make the assumption that 
mi << m2, the integrals Iq and I±, multiplied by 16n 2 , reduce to log(m|/mi) 
and m 2 //ii, respectively. We therefore find from (22) that 



With this choice of the behavior of the PV masses, the integrals involved 
in the structure functions (33) and (34) have no singularities (in m 2 ) worse 
than logarithmic. From (38) we see that if mi stays finite or diverges more 
slowly than m 2 /logm 2 in the large-m 2 limit, g will go to zero so fast that 
the structure functions must vanish, and we will have, in that sense, a trivial 
theory. A further examination of the particular choice m\ ~ m 2 / logm 2 shows 
that even in that case the structure functions go to zero as m 2 goes to infinity. 
The only choice for the behavior of m\ as a function of m 2 which leads to finite, 
nonzero structure functions is mi ~ m 2 . We therefore define r = m\jm 2 and 
hold r fixed as /i2 = m 2 ^ 00. The fractional amplitude ( for the single-PV- 
fermion state, given in (23), becomes equal to r. The coupling is then driven 
to a fixed value 



Thus mi (and r) must be negative. The structure functions become 




(m 2 - mi)[m 2 =F M log(mf /m|)] ' 



(M^mi)(M^m 2 ) 



(38) 




2 

r y 



1 + (r 2 - 1) y 

r ,(y) _ r v ^°s(r 2 y) 

y r 2 y _ I 
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r 2 y r 2 y r log(l — y + r 2 y) 



log(- + , , 

y 1 — y + r z y r A — 1 

ry l-y + r 2 y, 

log(- 



l + (r 2 -2)y &v y 



(40) 



f B .(y) - *(1 - r)(l - r 2 ) ( 4 y - (1 + r ^!° g(r2) + y(1+ 1 ^° g(y) 

| y (1 + r 2 y) log(r 2 y) y (1 + (2r 2 - 1) y) ^ r 2 y 

r 2 y — 1 1 — y 1 — y + r 2 y 

(2 + (r 2 -l) y) log(l -y + r 2 y) 
r 2 - 1 

+ l^+r^y (41) 

1 + (r 2 - 2) y y / 

The nonorthogonal projection of the wave function ensures that these distri- 
butions are positive definite. The reciprocal of the factor z\ is determined by 
the normalization condition (14) to be 



1 1 a 2 r 

_ = (1 - r) 2 + — -§— 4(27 - 108r + 307r 2 ) - 3tt 2 (6 - 24r + 49r 2 ) 



z\ v ' 14416tt 2 



-24(3 - 12r + 20r 2 ) ln(r 2 ) (42) 



to second order in r. We thus have a one-parameter family of theories labeled 
by r. While the PV masses have been taken to infinity, they have not been 
made infinitely large compared to the bare fermion mass, which has been 
taken to minus infinity. The value of g is finite in this limit. We probably 
should not, even naively, think that all the effects of the negatively normed 
states have been removed from the full solution. To control such effects, we 
should consider values of r which are small in absolute value. Notice that g is 
then restricted to small values. 

The results of the exact solution are very different from perturbation theory. 
In first-order perturbation theory, M is equal to mi, and there is no nonlinear 
eigenvalue equation and thus no restriction of the value of g. Indeed, since the 
physical mass, M, is fixed and equal to mi, we could not send mi to minus 
infinity as we did above. We also note that the discrete chiral symmetry is not 
restored in the large PV-mass limit. We cannot take mi to be zero (without 
obtaining a trivial theory). We can take the physical mass, M, to be zero, but 
that point does not occur at mi = 0. 

Plots of the structure functions for r = mi/m 2 = —0.01 are given in Fig. 3. We 
should remark on the behavior of the structure functions at the end points. For 
very large values of the PV masses the functions are given essentially exactly 
by Eqs. (40) and (41) for all points except very close to y = in the case of 



15 



The exact structure functions are zero at y = for all values of the PV 
masses; yet (40) yields a nonzero value at that point. Thus the convergence 
to the limiting forms is nonuniform. For that reason, any quantity sensitive to 
the endpoint behavior, such as the expectation value of the parton light-cone 
kinetic energy, should be calculated for finite values of the PV masses and 
then taken to the infinite-mass limit. 
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Fig. 3. Structure functions fs+ (solid) and (dashed) for the equal-PV-mass 
case, with r = m-\_/m,2 = —0.01, from the forms given in Eqs. (40) and (41) of the 
text. 



3.2 Unequal Pauli-Villars Masses 



Having obtained the results discussed in the last subsection, one can ask 
whether there is any way to get results more like perturbation theory. As 
it turns out, there is: to do so we must take the limit of large PV masses in 
such a way that the PV fermion mass grows much faster than the PV boson 
mass. 4 

If we take the mass m 2 to infinity, the integrals Iq and I\ reduce to 

167T 2 / ~ - log^/^O , 167T 2 /! ~ -2— log^/Vi) . (43) 

Hi 



The fractional amplitude ( = (m x =f M)/m 2 goes to zero. When we take 

4 We could let the boson mass grow as fast as logrn^, but it is also allowed, and 
is simpler, to first take ni2 to infinity at finite ^2 (that limit turns out to be finite) 
then take ^2 to infinity. 
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m 2 — > oo and then take /i2 large, the eigenvalue equation (22) becomes 

9 2 _ C 
16tt 2 log(/i 2 //ii) ' 



(44) 



where 



C=( mi ™ ). (45) 



In this limit the structure functions reduce to 



fB+{y) 



f B -(y) 



g 2 4 



[l±m 1 /ii 1 -y) 2 yiA 



16n 2 —y)+ m\y 
g 2 z\ 



16tt 2 



y { log 



M 2 y{\ - y) ' 
(i - 2/ VI 



A^K 1 - y) + m \y - M 2 y(l - y) 



- 2 



(46) 
(47) 



3.2.1 mi finite 

Looking at these relations we see that if g 2 ~ 1/ log/Z2, /s- will be finite and 
nonzero while J'b+ will be zero in the limit of large PV mass. If g 2 remains 
finite, /b+ will be finite and nonzero while will diverge, which is untenable. 
There are two choices for the behavior of mi which will give us the desired 
behavior for g 2 and finite non-zero values for fs-iy). One way is to choose m\ 
to be finite and choose its value and the signs in (45) such that the constant 
C is any real number we wish. In that case, fs- is given by 

f B -{y) = 2Cz\y. (48) 



From (26) we find that 



Thus there is a finite probability that the state consists of a single physical 
fermion. The larger the value of C the smaller is that probability and the 
larger is the probability that the state contains two particles. We note, as in 
the case of equal PV masses, that the discrete chiral symmetry is not restored 
in the sense that if we take either M or mi equal to zero, the other is not 
specified and disappears entirely from the problem; the value of C is fixed at 
either 1 or 1/2. 
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3.2.2 mi proportional to M 



The other possibility for the behavior of m\ is to choose m\ ~ ±f with the 
appropriate choice of sign in (45). For illustration we take the lower sign and 
parameterize 

mi = T + 2clogW^)' (50) 



with c a constant. With these choices the bare coupling constant goes to a 
finite value given by 

167T2- 2 • (51) 



From this we see that c should be positive. Notice that this choice is much more 
like perturbation theory: instead of mi = M, as in first-order perturbation 
theory, we have (in the limit) mi = : y, and the coupling constant can be any 
finite number. In this case we find for large \ii that z\ is given by 

2 

Z ' = 3Mclog(^M ■ (52) 



There is zero probability that the system is in the state of one physical fermion, 
and the entire wave function is in the two-particle sector. Due to the behavior 
of z\ we find again that, in the infinite-/^ limit, Jb+ is zero while 



/fl-(y) = 2y. (53) 

The outcome for the discrete chiral symmetry in this case is not so clear. 
The fact that mi is proportional to M (in the limit) suggests that it may be 
restored. On the other hand if M is zero we would encounter undefined expres- 
sions in the above derivation. However, we can perform the entire calculation 
with M set equal to zero from the start, and we find that if we take 

m >=2£^- < 64) 



we obtain the structure function (53) and zero for f B+ (y); this last result is 
in agreement with perturbation theory. So in that sense, the discrete chiral 
symmetry is restored in the large-PV-mass limit. 

We should repeat the comment of the previous section regarding the behavior 
of the structure functions at the endpoints. For finite values of the PV masses, 
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the structure functions vanish at y — 1, but there is a nonuniform convergence. 
For very large values of the PV masses the structure function is closely pro- 
portional to y for all values of y except very near 1 where it falls precipitously 
to zero. In the limit of large PV masses the function converges to something 
proportional to y for every point except y — 1, where it is always zero. For 
that reason any quantity which is sensitive to the endpoint behavior (such as 
the kinetic energy of the fermion or (:</> 2 (0):)) should be calculated for finite 
values of the PV masses then the limit taken. If that exercise is performed for 
(:0 2 (O):), we find that this quantity does diverge. 



4 On Not Taking the Limit 



Up to now we have taken the limit of the PV masses going to infinity. Here 
we wish to further consider the comparison of our results with perturbation 
theory. We believe that this comparison suggests that we should not take 
that limit and furthermore indicates why we should not do so. These same 
considerations will suggest a way to decide how large we should take the PV 
masses. 

Let us fix our attention on the choices made in Sec. 3.2.2 for taking the limit 
of large PV masses and fixing mi, which gave results most like perturbation 
theory. The structure function was zero in that case. That does not 

happen in perturbation theory. Since our wave function is identical and even 
the parameters are almost the same (differing only in that mi = M/2), how 
can we get something so different from perturbation theory? The reason we 
obtained zero for fs+ is that the renormalization constant Z\ went to zero. If 
we look at the form of the function for finite PV mass, it is 

f B ( y ) = ^ 2 [finite quantity] 

1 + g 2 [finite quantity] + g 2 [finite quantity] log/x 2 



The denominator represents z{ 2 . Now in perturbation theory, since the nu- 
merator is already of order g 2 , only the 1 in the denominator is used and the 
result is nonzero. Indeed, suppose we calculate some quantity which is finite 
to this order such as the anomalous magnetic moment. Again we would get a 
result of the form 

= g 2 [finite quantity] 

1 + g 2 [finite quantity] + g 2 [finite quantity] log /i 2 



If we use the methods of the previous section this quantity would again be zero. 
In perturbation theory that would not happen, again because the divergent 



19 



term in the denominator would not be used with this numerator. Now the 
divergent term in the denominator would be used in a calculation to order 
g A ; but then there would be an order g 4 term in the numerator which would 
cancel the divergence of the term from the denominator and give a finite result. 
That is the way perturbation theory works. The point is this: we will have an 
accurate calculation only to the extent that the projection of the wave function 
onto the excluded Fock states is small. We know from past calculations [10] 
that this projection can be very small, sometimes even for the severe truncation 
we are considering here, but those results were for finite values of the PV 
masses. There will be divergences in the excluded Fock sectors, and we must 
anticipate that for sufficiently large values of the PV masses the projection of 
the wave function onto those sectors will not be small. 

There are two types of error associated with having finite values of the PV 
masses: for PV masses too small we will have too much of the negative norm 
states in the system. We anticipate that such errors are approximately the 
larger of — or where mp is the smallest PV mass. The other type of error 

° mp rap r J 1 

is a large projection of the wave function onto the excluded Fock sectors. That 
error should be approximately 

\ +physl^+phys/ 
($ + p h ys|$+phys) 



where 3>+ p hys is the projection of the wave function onto the lowest excluded 
Fock sector. 5 The higher Fock wave function $+ p h ys can be estimated using 
perturbation theory, perturbing about $ + with the projection of P~ onto the 
excluded sectors being chosen as the perturbing operator. The first type of 
error, from negative-metric Fock states, decreases with increasing PV mass; 
the second type of error, the truncation error, will usually increase with in- 
creasing PV mass. Ideally we should choose the values of the PV masses to 
be the values where the two types of error are equal. The strategy for treating 
the nonperturbative system is to include more and more of the representa- 
tion space in our calculation and to increase the value of the PV masses until 
the desired accuracy is achieved. How much of the space will be required will 
depend on the problem. 

In later work we will attempt to make these comments quantitative by es- 
timating the optimum values for the PV masses. Here we will illustrate the 
effects of not taking the limit, for the trajectory in which the PV masses are 
taken to large values as in Sec. 3.2.2. For infinite m-i but finite fi2 (= 100/xi), 
and with the eigenvalue equation solved by Eq. (50) and Eq. (51), the struc- 
ture functions, fs+ and are plotted in Fig. 4. These are to be compared 

some rule other than particle number is used to truncate the space, ^^p^yg is 
the projection onto the "next" set of vectors. 
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with the linear function (Eq. (53)) for fs-, and zero for fs+, which result if 
the limit of infinite PV masses is taken. The structure functions for the finite 
value of the mass M resemble what one expects in a bound state. 
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Fig. 4. Structure functions (a) and (b) /b+ for the unequal-PV-mass case 

where 1712 —>■ oo, ^2 = 100/xi, and M = jjL\. The bare fermion mass mi is specified 
by Eq. (50) of the text, and the coupling g 2 by Eq. (51), with c = 1. Notice that 
the two plots have different vertical scales. 



5 Discussion 



In this paper we have studied the regulation of Yukawa theory by the use of 
Pauli-Villars fields in such a way that the interaction is written as a product 
of zero-norm fields. Paston and Franke have shown that this regulation proce- 
dure gives perturbative equivalence with Feynman methods [12]. The theory is 
covariant and presumably finite, and there are no gauge symmetries to protect. 
Therefore, if we could solve such a theory exactly and take the limit of the PV 
masses going to infinity, the result would be the best one could do to give a 
meaning to the theory. In this paper we have done our calculations in a severe 
truncation of the representation space. Such a truncation violates covariance, 
but if the contribution of excluded Fock sectors is sufficiently small, the con- 
sequence of the truncation is more a question of accuracy than of preserving 
symmetries: if we are close to the hypothetical solution mentioned above, it 
does not matter if this (small) error violates symmetries. 

The reason we have used such a severe truncation is that it allows us to find 
solutions and take limits in closed form, and thus our interpretation of the 
results is not confounded by questions of inaccuracies introduced by numerical 
solutions. A significant feature of the calculations, which came as a surprise 
to us, is that the results depend strongly on the way in which the two PV 
masses are allowed to approach infinity. It is not true that any two different 
trajectories will give different results but rather that there are families of 
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trajectories which give the same results. For instance, any trajectory on which 
the two PV masses are proportional to each other, with a fixed constant of 
proportionality, give the same result as taking the limit with the two masses 
set equal to each other. Similarly, any path on which /i 2 is logarithmically small 
compared with m 2 will give the same result as taking the limit m 2 — > oo first, 
then taking the limit \x<i — > oo. One possibility is that all these trajectories 
represent different phases of the theory. Another possibility is that the effect 
is an artifact of the truncation, and, if we include more and more of the 
representation space in the calculations, the results of the various ways of 
taking the limit will approach each other. Another possibility that we have 
considered is that some of the ways of taking the limit are wrong and that 
some principle which we have not yet discerned will determine the correct way 
to take the limit. We hope to report further studies on this question in the 
future. 

In the calculations we have given special consideration to the discrete chiral 
symmetry that is formally present in the unregulated Lagrangian. Writing 
the interaction as a product of zero-norm fields breaks the chiral symmetry 
explicitly (unless the mass of the PV fermion is taken to be zero), and we have 
been careful to notice whether or not it is restored in the large-PV-mass limit, 
at least in the sense that the physical mass of the fermion is proportional to the 
bare mass. We find that for some ways of taking the limit the discrete chiral 
symmetry is restored and for some ways it is not. We do not know whether 
this consideration can provide a valid way of choosing one limiting procedure 
over another. This question is important because, not only is chiral symmetry 
of interest in itself, but the way it is broken by the regulation procedure is 
analogous to the way gauge symmetry is broken by writing the interactions of 
gauge theories as products of zero-norm fields. 

We have argued that our results suggest that if calculations are done in a 
truncated representation space, it may not be correct to take the limit of the 
PV masses going all the way to infinity. It is easy to understand the reason 
why: if we are to have accurate calculations, most of the support of the wave 
functions in which we are interested must lie in the part of the space we retain. 
We know from past studies that the projection of the low-lying states onto the 
higher Fock sectors often falls off very rapidly in the light-cone representation, 
but those results were for finite values of the regulators. At infinite values of 
the regulators, the eigenvectors are not expected to exist at all, and we must 
expect that as the regulators are removed the projection of the wave functions 
onto any allowed sectors will become large. Thus it will be necessary to keep 
the PV masses finite when one truncates the representation space. If there are 
values of the PV masses sufficiently large to remove most of the bad effects of 
the negative-norm states on the eigenvectors in which we are interested, but 
small enough to make small the projection of these eigenvectors onto sectors 
we cannot manage to keep, then we can do a useful calculation; otherwise not. 
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We are currently performing studies to try to make these remarks quantitative. 
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